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Abstract 

The theoretical formulation of x-ray resonant magnetic scattering from rough surfaces and in- 
terfaces is given for the diffuse (off-specular) scattering, and general expressions are derived in 
both the Born approximation (BA) and the distorted-wave Born approximation (DWBA) for both 
single and multiple interfaces. We also give in the BA the expression for off-specular magnetic 
scattering from magnetic domains. For this purpose, structural and magnetic interfaces are de- 
fined in terms of roughness parameters related to their height-height correlation functions and the 
correlations between them. The results are generalized to the case of multiple interfaces, as in the 
case of thin films or multilayers. Theoretical calculations for each of the cases are illustrated as 
numerical examples and compared with experimental data of mangetic diffuse scattering from a 
Gd/Fe multilayer. 



I. INTRODUCTION 



In the preceding paper,— we have developed the dynamical theory for x-ray resonant 
magnetic specular reflectivity using the self-consistent method in the distorted-wave Born 
approximation (DWBA). It is important to bear in mind, however, that specular reflec- 
tivity, which measures the density profile normal to the surface averaged over the in-plane 
directions, can yield only information corresponding to the root-mean-square roughness (or 
equivalently, the average width) of the interfaces. A more complete description of the mor- 
phology of the roughness can only be obtained from off-specular or diffuse scattering studies. 
The first such studies of resonant x-ray magnetic diffuse scattering studies were carried out 
recently by MacKay el al£ From these measurements, quantities, such as the in-plane cor- 
relation length of the roughness, the interlayer roughness correlations and the roughness 
exponent, can be deduced. These quantities are of considerable importance. In the case 
of magnetic films, Freeland et a/.— found little correlation of the variations in the magnetic 
coercive force H c for a variety of samples with the average roughness (chemical or mag- 
netic) but a systematic dependence on the roughness correlation length. In the case of giant 
magnetoresistance (GMR) films, contradictory results have been found in studying how the 
magnitude of the GMR effect depended on the chemical roughness alone, and it is likely that 
the effect depends on a more detailed set of parameters related to the magnetic (as well as 
possibly the chemical) roughness.^ Barnas and Palasantzas^ have carried out calculations 
of the manner in which self-affme roughness at an interface affects electronic transport. 

Methods were developed earlier to calculate analytically the offspecular components of the 
charge scattering of x-rays by rough surfaces and interfaces using the Born approximation 
(BA) and the distorted- wave Born approximation (DWBA).— We present here the gener- 
alization of these methods to the case of resonant magnetic x-ray scattering from surfaces 
or interfaces of ferromagnetic materials possessing both structural and magnetic roughness. 
For this purpose, in the preceding papepi we have represented the deviations from a smooth 
magnetic interface in terms of a "rough" magnetic interface, distinct from the structural 
interface (but possibly correlated strongly with it), with its own self-affme roughness pa- 
rameters and parameters representing the correlation of the structural with the magnetic 
roughness height fluctuations. Components of the magnetization at the interface, which are 
disordered on much shorter length scales, are ignored in this treatment, as they will scatter 
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at much lager q-values than those of interest here. The BA results have been previously 
presented in an earlier publication^ and already applied in interpreting x-ray resonant mag- 
netic diffuse scattering measurements from magnetic multilayers.— We also note that the 
analogous case of off-specular neutron scattering by magnetic roughness was treated earlier 
by Sinhaifi' and more recent treatment has been given by Toperverg.— 

The plan of this paper is as follows. In Sec. II, we derive the expression for the magnetic 
scattering in the BA, which is presented here in detail for completeness, although a brief 
account has been published earlier.— In Sec. Ill, we discuss the magnetic diffuse scattering 
from magnetic domains in the BA. In Sec. IV, we present the derivation of the resonant 
magnetic diffuse scattering in the DWBA for a single magnetic interface and discuss numeri- 
cal results. Finally, in Sec. V, we discuss the extension of the formalism to the case of diffuse 
scattering from magnetic multilayers and present some numerical results with experimental 
data from a Gd/Fe multilayer, which were analyzed earlier in the BA.— 



II. RESONANT MAGNETIC X-RAY SCATTERING IN THE BORN APPROXI- 
MATION 



In Sec. Ill of paper I,— we discussed the dielectric susceptiblity Xa/3 of a resonant magnetic 
medium. The resonant scattering amplitude density can also be given by 

k 2 

= [-^oPo(r) + n m (r)A) 5 a/3 

- iBn m {v) ]T e a/?7 M 7 (r) + Cn m (r)M a (r)Mp(r), (2.1) 
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where a, (3 denote Cartesian components, A, B, C are the energy-dependent parameters 
defined in Eq. (2.3) of paper I,— and M a is the a-component of the unit magnetization vector 
of magnetic atoms. We may write an effective total electron number density function 



Peff(r) = Po(r) n m (r) 

?"0 



and write 



•7^/3 (r) = -r 



Peff (r)<5 a/3 + iBn m (r) e Q/?7 M 7 (r) - C , n m (r)M a (r)M /3 (r) 



(2.2) 



(2.3) 



where B = B/r , C = C/r . 



Let us consider the cross section for scattering for a photon from a state |kj, \i > to a 
state |k/, v >, where (kj, fi), (k/, v) represent, respectively, the wave vector and polarization 
state of the incident beam and those of the scattered beam. This is given by 

© k k =i^i<^mw>i*, (2.4) 

where < • • - |T| • •• > denotes the matrix element of the scattering from the interface (s). Let 
us define shown schematically in Fig. 1, where the z-axis is normal to the average 

plane of the interface. In the Born approximation (BA), the matrix element can be written 

as 

< k f ,i/\T\\£i,n>= 4vr^e* Q e M/3 / dre^F^), (2.5) 

where e u , are the photon polarization vectors corresponding to the final and incident 
photon state, respectively, and q = kf — kj. Using Eq. ()2.3|) . we obtain 



<k f ,u\T\ki, fx > = -47rr 



+ * E (q) - £ et a e^CM^ (q) 

a/37 a P 

where 



(2.6) 



Pefr(q) = J dre iq ' r p eff (r), 
M«(q) = J rfre- iq - r n m (r)M 7 (r), 

M$(q) = J dve-^ T n m {v)M a {v)Mp{r). (2.7) 

We now restrict ourselves to the simplified model discussed in Sec. II of the paper I,— 
of a single interface between a magnetic and a nonmagnetic medium with a chemical (or 
structural) interface defined by a height z c (x, y) and a magnetic interface defined by a height 
z m (x,y). In accordance with previous approximations valid for small q <C a -1 , where a is 
an interatomic spacing, we may assume p e ff(r) constant with the value p\ for z < z c (x,y) 
and the value p 2 for z > z c (x,y), and n m (r), M(r) constant for z < z m (x,y) and zero for 
z > z m (x,y). In this case, after carrying out the ^-integration, Eq. ()2.7|) simplifies to 

p cff (q) = / Pl ~ p2) f f dxdye~ ici r p e- iq * Zc{x ' y \ 
q z J J 

M «(q) = i 1 ^! [ [ dxdye' l ^P e ' iq ^ y \ 
7 Qz J J 

r (2), s _ , n m M a M p f f ■ - igzZmix , y) 



Af$(q) = i - - - p / / dxdye-^Pe-^^y), (2.8) 
Qz 



where qn and p are in-plane components of q and r, respectively. Let us denote the heights 
z c (x,y), z m (x,y) collectively as Zi(x,y) [i = c,m] and define 



G c = (pi-p 2 )(e* e 



B(et x e„) • M + iC{e* v ■ M)(e M • M) . (2.9) 



Then from Eqs. ()2.4|) and ()2.6|) we obtain 

(s) =? E GiG*Sij, (2.10) 

where 

S iJ = J J J J dxdydx'dy'e-^p-p')e- iq ^ x ' y) - z > (x '> y,)) . (2.11) 

Except in cases where the incident x-ray beam is coherent over the sample, we may 
introduce the usual statistical configurational averages to evaluate Sy to obtain 

Sij = Ae- iq ^-~ z ^ f f dXdYe-^/e-^ 6 *^- 5 ^ ' ^), (2.12) 



where A is the illuminated surface area, R [= (X,Y) = (x — x',y — y')] measures the in- 
plane separation of two points on the appropriate interfaces, and 8z c (R), 5z m (R) are the 
height deviations from the average heights z c , z m of the chemical (structural) and magnetic 
interfaces, respectively. We allow for the possibility of a finite separation A between the 
structural and magnetic interfaces, as shown in Fig. 1. If we make the customary Gaussian 
approximation for the height fluctuations Sz c (x,y), Sz m (x,y), respectively, we obtain 



where 



S i3 = Ae-^-^ J J rfXrfFe- iq ii- R e-5^([^( x ^)-^(°' )] 2 ) 
= z-WM+^Ae-^-^ J J dXdYe-^e^W 

= e-^^ + ^ ) e- iq ^-^hn 2 A5(q x )5(q y ) + SL = c,m) (2.13) 



S' {j = e'hl^+^e-ii^-^A \ \ dXdYe- 1 ^ e & c vW - 1 , (2.14) 



C CC (R) = < 6z c (0)6z c (R) >, 
C mm (R) = < 5z m (0)5z m {R) >, 

C cm (R) = < 6z c (0)6z m (K) >, (2.15) 



and a c , a m are the root-mean-squared chemical (structural) and magnetic roughnesses, re- 
spectively. Equation (j2.13|) contains expressions for both the specular scattering (arising 
from the 5- function containing terms in S" _, ) and the diffuse (or off-specular) scattering, 
which include both charge and resonant magnetic scattering. 
Collecting the specular terms, we obtain 

fda\ spec _ 47rUrg 



- ) = —^S(q x )S(q y )Q, (2.16) 



where 



Q= J2 e-i« , 'tf +a j ) e- iq *te-*t)GiG* j , (2.17) 

i,j=c,m 

from which the specular reflectivity can be immediately expressed as Rn^v = — t^Q- The 
diffuse scattering may be expressed as 

{%) =? £ GiG '* s «' (2 ' 18) 

V / ^— >i/ i,j=c,m 

The explicit expressions for the specular reflectivity and diffuse scattering from a rough 
surface in the BA for specific directions of the magnetizations and the photon polarizations 
are given in Appendix A. 

Let us now consider a multilayer with N interfaces (N — 1 layers). Each layer can be 
characterized by its effective total electron number density p n for each layer, and by resonant 
magnetic scattering amplitude densities n m ^ n B n , n m ^ n C n and magnetization vector M n for 
resonant magnetic layers. The n-th interface lies between the n-th and (n + l)-th layers, 
and its average height is denoted by z n , as shown in Fig. El 

The specular reflectivity from a multilayer with N interfaces can be then expressed as 

iW = ^^£ E e-^^e^^-'^Gt,^, (2.19) 

^2 n,n' i,j=c,m 



where 



G c<n = (Pn+1 ~ Pn) (K ■ e u 

G 



-in r , 



B n+1 (e* u x e M ) • M n+1 + iC n+1 (e* u ■ M n+1 )(e M • M„ +1 ; 
B n {e* v x e„) • M n + iC n (e* u ■ M n )(e„ • M n )l. (2.20) 
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The diffuse scattering from a multilayer may be also expressed as 

/dcr\ difrusc r 2 N 



jqJ " ? ° 2 E Gi,nG* j n/ S* jnn ,, (2.21) 



where 



„2 



f2.22) 



and the cross-correlation functions between the n-th and n'-th interfaces are defined by 

C cc , nn /(R) = (8z Cin (0)5z Ctn >(R)) , 

C cm>nn ,(R) = (6z Cjn (0)5z mtn i(R)) . (2.23) 

Here 5z CtU , 5z m ^ n are the height deviations from the average heights of the n-th structural 
and magnetic interfaces, respectively. We also allow for the possibility of a finite separation 
A n between the n-th structural and magnetic interfaecs, as shown in Fig. 2. Explicit forms 
for these cross-correlation functions have been given by several authors , 8 i 9 i 12 i 13 and may be 
substituted in Eq. ()2.22|) to yield the diffuse scattering cross sections. Numerical examples 
for the reflectivity and magnetic diffuse scattering for single surfaces and multilayers have 
been given in previous publications.^ We defer showing these here until Sec. IV, where 
we present them together with the results from the distorted-wave Born approximation 
(DWBA). Explicit forms for specular reflectivity and diffuse scattering from multilayers for 
specific directions of magnetizations and the photon polarizations are also given in Appendix 
A. 

The Born approximation also allows one to include explicitly the effects of a graded rather 
than sharp (but rough) magnetic interfaces, i.e., where the magnetization in the resonant 
medium is not uniform but decays towards the interface. Let us suppose that in the vicinity 
of a particular interface, the magnetization component M a (r) decayed with an envelope 
function $>(x,y, z) (we assume this is independent of the component a). The height of the 
magnetic interface z m (x,y) can then be defined at the lateral position (x,y) as the position 
of the point of inflection of $, i.e., where t^<&(x, y, z) = 0. 

If we assume that the function y, z — z m (x,y)) is independent of (x, y) and can be 
written as <&(A#), then we can define a form factor ip(q z ) by 

cp(q z ) = J d(Az)^^-e-^ Az . (2.24) 
7 



Then the effect of the graded interface can be introduced by simply multiplying S cm in Eq. 
(I2.12J1 by <f(q z ) and S mm by \{p(q z )\ 2 . The generalization to the case of multiple interfaces 
in Eq. (|2.22|) is obvious. 

We should mention, however, that in this case a m in Eqs. (J2.12j) and (|2.22j) represents 
purely rough interfacial width rather than the total interfacial width at the magnetic in- 
terface. The latter includes the effects of both graded interface due to interdiffusion and 
interfacial roughness, and specular reflectivity discussed in the preceding paper 1 provides 
only this total interfacial width. On the other hand, since the correlation functions in Eqs. 
()2.15|) and ()2.23j) contains only pure roughness diffuse scattering allows one to distinguish 
pure interfacial roughness from the graded interface.— 



III. MAGNETIC DIFFUSE SCATTERING FROM MAGNETIC DOMAINS IN 
THE BA 

For an unmagnetized or patially magnetized film, there will exist magnetic domains which 
can give rise to off-specular (diffuse) magnetic scattering even in the absence of magnetic 
roughness. There are two cases to consider: 

(I) If the typical lateral size of the domains is larger than the lateral coherence length 
of the x-ray beam on the surface of the film, then the scattering will be the sum of the 
scattering from the magnetized regions. If there is no net magnetization, the terms linear 
in the magnetization which appear in the interference between the charge and magnetic 
scattering in Eqs. (j2.18J) and (jA9|) will cancel out. In particular, if we neglect the other 
terms [without C in Eq. (jA9|) ] there will be no magnetic contribution to (f§|) + — (^) 
Domain scattering will not manifest itself in the off-specular scattering as the length scale 
invloved will be too large for the qn at which it would occur to be resoved from the width 
of the specular reflection. 

(II) If the lateral size of the domains is smaller than the lateral coherence length of the 
x-rays, domain scattering will manifest itself in the off-specular scattering. (This will also 
be true of magnetic clusters, dots or other laterally inhomogeneous magnetic structures in 
the film). Let us consider a simplified model where the domains are parallel or antiparallel 
to the average direction of magnetization. Ths can be expressed by a function p(x, y) which 
takes the values of +1 for the domain magnetization being parallel to the average direction 
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and —1 for the domain magnetization being antiparallel. 

Then if we assume that the function p(x, y) is uncorrelated with structural features at 
the interface such as the structural rouhness, we can define a statistical two-dimensional 
domain correlation function 

ld {X, Y) = (p{x, y)p{x + X,y + Y)), (3.1) 

which is analogous to the Debye correlation function for a three-dimensional porous medium. 
This is related to the probability that a vector (X, Y) on the surface has one end in one 
domain and the other end in a similarly oriented domain. Then the expression for Sy in Eq. 
()2.12j) with i,j = m must be modified to 

S mm = A[[ rfXrfFe- tq ii- R (e^ (<5z '" (R) ^ m(0 ' 0)) )7 d (X,F), (3.2) 



while S cc , S cm are unmodified. In particular, if we neglect magnetic roughness and consider 
"only" domain effects, 

S mm = A [ [ dXdYe-^ ld (X, Y), (3.3) 



which will give rise to magnetic domain scattering from an otherwise smooth surface. 

A common approximation for random domains and sharp walls Jd{X, Y) = e~ R ^ a , where 
a is the average domain size, yields 

a 2 

(1 + q^a 2 y/ 2 

From Eq. ()2.1())) we see that such domains will give magnetic off-specular scattering which 
will behave asymptotically as q^ 3 which is the two-dimensional analogue of Porod scattering 
from random smooth interfaces. 

IV. DIFFUSE SCATTERING FROM A SINGLE MAGNETIC INTERFACE IN 
THE DWBA 

The diffuse scattering in the distorted-wave Born approximation (DWBA) can be given, 
from Eq. (4.13) in the preceding paper,— by 

( da\ 1 



- 



dVt ) 167T 2 

% / dittusc 



T fl \ 2 ) - (T fi " 2 



(4.1) 



where T? % =< k/, z/|T|kj, // > is the scattering matrix element with the vector fields |kj, /i > 
and |ky, v >, which were defined in Eqs. (4.4) and (4.10) of paper I.— As shown in Eq. 
(4.14) of paper I,— < kj, z/|T|kj, /i > can be approximated in the DWBA by the sum of three 
matrix elements involving x \ A c , an d A m , which represent an ideal system with a smooth 
interface and perturbation on x due to structural (chemical) and magnetic roughnesses, 
respectively, as defined in Sec. IV of paper I — Bearing in mind that the first term involving 
X^ vanishes for diffuse scattering, the remaining terms involving only perturbations A c and 
A m contribute to diffuse scattering, and their matrix elements are evaluated with the vector 
state |kj,/i > [or E(kj,/i) in Eq. (4.4) of paper Ji] and its time-reverse state | — kj, v > 
[or E T (— ky,i/) in Eq. (4.10) of paper ll]. Assuming for both E(kj 5/ u) and E T (— ky,z/) the 
functional forms for z < analytically continued to z > 0, in the spirit of Ref. |7j, gives 

<k f ,v\T\Kn> = E^ 0) (-k/)^(k,^ iK 



x 



Qtz(jj') 

(Xi -Xo)E e ,W« / / dxdye-^'^^e 



jj 

igtz(jj')Sz c (x,y) „-iqii p 



+ J2 e jaX ( Se ff , I I dxdye-^^ 5z ^e-^P 

af3 



(4.2) 



where 



q i (j/) = kUj)-k*(j / ), (4-3) 



and Xo,i an d Xal are dielectric susceptibilities defined in Eq. (4.2) of paper li, and k*(j') 
and kj(j) are the transmitted wave vectors for the incident and scattered waves, as defined 
in Fig. 2 of paper I.— Here (— k/) and T^(kj) in Eq. ()4.2)1 are 2x2 matrices denoting 
the transmission coefficients for (— kj) and kj waves in Fig. 1, respectively, from the average 
smooth interface. Their explicit forms for small angles (Of -C 1 for the incidence angle 9i) 
and M || x were given in Appendix A of paper I.— We note that because of the continuity 
of the wavefields [see Eqs. (4.5) and (4.11) of paper I 1 ] 

(q*(j/))|| =Q|| for all (4.4) 

Substituting in Eq. (|4.1jl and carrying out the statistical average over the interface, we 
obtain 



- , \ diffuse , a 

da \ fen 



E T$\-*t)T${-±t)TftM^ (4.5) 
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where 



+ { (Xi - Xo)* (I] e iQ e*, a ) (J2 elyX^ek's) U cm + c.c. } 



(4.6) 



and 



Uw = J J J J dxdydx'dy' 



In Eq. ()4.6j) c.c. refers to the complex conjugate of the preceding term, and in Eq. ()4.7j) we 
have again made the customary Gaussian approximation for 5zi(x,y) [I = c,m] in carrying 
out the average. The notations are the same as in Eqs. (|2.1H|) - (l2.15j) . 

For circularly polarized incident x-rays with e±(kj) = [e cr (kj)±ie 7r (kj)]/ y/2, the scattering 
intensities without polarization analysis for the outgoing beam can be evaluated as 

© ± = ih\ i:j<*,MTK°>±i<* f MT\K*>\\ (4.8) 



+ 



X 



T^(k i )7]g(k i )H( & (7jO ) &(A ; AO)}. (4.9) 



We shall now illustrate numerical examples calculated again for a Gd surface with various 
structural and magnetic roughnesses and magnetizations along the x-axis. Here, we have 
used the height-height correlation functions introduced by Sinha et al. 7 for self-affme fractal 
interfaces, i.e., 

C cc (R) = a c V^) 2,lc , C mm (R) = ale~& hm , (4.10) 

where £ c and £ m are the lateral correlation lengths of the roughnesses at the structural and 
magnetic interfaces, respectively, and h is the roughness exponent describing how jagged 
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the interface is. For the cross-correlation function between the structural and magnetic 
interfaces separated by a magnetically dead layer, as shown in Fig. 1, the Schlomka et a/.— 
expression was used: 

C cm (R) = — — (e +e y irn> j e e J-.«», (4.11) 

where £_i_ jCrn is the vertical correlation length between the structural and magnetic interfacial 
roughnesses separated spatially by A. 

Figure El shows the calculations of x-ray resonant magnetic diffuse scattering intensities 
at the Gd L 3 -edge from Gd surfaces with different roughness parameters. In Fig. Efa) and 
(b), the structural interface has a roughness of a c = 5A and £ c = 800A, and the magnetic 
one has a m = 3A and £ m = 1500A. On the other hand, in Fig. |3fc) and (d) the values of 
roughness parameters for the structural and magnetic interfaces were reversed, and those 
interfaces are separated by a 20-A-thick magnetically dead layer, as shown in Figs. 3(g)- 
(i) of paper I.— The roughness exponent h = 0.8 was used for all structural and magnetic 
interfaces. Longitudinal diffuse scattering intensities in Figs. Efa) and (c) show similar 
features as specular reflectivites in Fig. 3(a) and (g) of paper Ji and do not show clearly the 
effect of different lateral correlation lengths. Instead, since longitudinal diffuse scattering 
intensity is sensitive to the vertical correlation length between spatially separated interfaces, 
for the dead-layer sample we performed calculations with and without vertical correlations 
between the structural and magnetic roughnesses separated by the magnetically dead layer. 
These vetical correlations are related mainly to the interference term, (I + — /_), as shown 
in Fig. Efc), where circles represent a complete vertical correlation = 00). 

In order to show the effect of the lateral correlation lengths, we have also calculated 
transverse (or rocking curve) diffuse scattering intensities, where the widths of the diffuse 
parts depend primarily on the lateral correlation lengths £. Figure Efb) and (d) show 
the transverse diffuse scattering intensities at q z = 0.2242A -1 normalized to unity by the 
maximum diffuse scattering intensities to clarify the effect of lateral correlation lengths. 
From two opposite cases, Fig. |3fb) and (d), with reversed values for the structural and 
magnetic lateral correlation lengths, we find that the width of the diffuse part of each 
scattering channel can clearly give a direct estimation of the corresponding lateral correlation 
length, as discussed above, i.e., a — > a scattering (solid lines) vs £ c and o — > n one (dashed 
lines) vs £ m . The widths for (7 + — /_) correspond to the effective lateral correlation length 
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£ cm between the structural and magnetic interfaces, which can be defined from Eq. (j4.11j) 
by 

(4.12) 

and then should be between £ c and £ m . This can be clearly seen in Fig. [2Kb) and (d). 

V. DIFFUSE SCATTERING FROM MULTIPLE MAGNETIC INTERFACES IN 
THE DWBA 

For a multilayer with multiple interfaces, following the representation for a single surface 
taken in Sec. IV of paper I,— each layer can be characterized by its dielectric susceptibility 
tensor Xa/3,n for the n-th layer, which can be Xa/3,n = Xn&a(3 for nonmagnetic (isotropic) 
layers and Xaf3,n = Xn^ap +Xaln f° r magnetic (anisotropic) layers, as defined in Eq. (3.5) of 
paper I.— The solution for the electric fields inside the n-th layer in the case of the "smooth" 
interfaces can be given by 

E n (k,, u) = £TiS(k 4 , u)e,- n e*^> r + E R^k,, u)% n e^>>, (5.1) 

3 3 

where the amplitudes T^°j(kj,/i) and R^j(k i} /i) are the vectors (T^l,T^l) and (R^^R^) 
representing two transmitted and two specularly reflected waves with wave vectors kj n (j) 
and k. r in (j) in the n-th layer, respectively, excited by an incident wave in vacuum with wave 
vector kj and polarization /i. These amplitudes can be obtained from recursive 2x2 matrix 
formalism developed by Stepanov and Sinha 15 and their explicit expressions were given in 
Appendix E of paper I.— The index j represents two components of each field amplitude and 
defines a-, 7r-component for nonmagnetic layers and (1)-, (2)-component for magnetic ones, 
as shown in Appendix A. Similarly to Eq. (J5.1)) . the time-reversed waves in the n-th layer 
for an incident wave with vector (— kf) and polarization v can be also given by 

El(-k /; = £ Tif(-k f , /^,., ( ' k - ^ - + E <^j,y K ' :j - r - (5.2) 

3 3 

For the structurally and magnetically rough interfaces in multilayers, we may write more 
generally 

N 

XafiW = E (xSLtt + A^ tB (r) + A^(r)) , (5.3) 

n 
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exp 




cxp 




+ exp 




where 



Xaln( r ) = Xaf3,n, for Z n < Z < Z n _l, 

= elsewhere, (5.4) 



A a/3,n( r ) = (Xn+i ~ Xn)8<xf3, for z n < z < z n + 5z c>n (x,y) if 5z CtU (x,y) > 0, 
= -(Xn+i ~ Xn)S a 0, for z n + 5z c<n (x, y) < z < z n if 5z c , n (x, y) < 0, 
= elsewhere, (5.5) 

and 

A a/3,n( r ) = (Xaln+1 ~ Xa},n)» for Z n < Z < Z n + 8z m>n {x,y) if 8z m>n (x,y) > 0, 
= -(Xaln+1 ~ Xahn)> ?n + 5z m , n (x, y) < Z < Z n if 5z m , n (x, y) < 0, 

= elsewhere, (5.6) 

5z C:n (x,y) and 5z mtn (x,y) denote the height deviations of the n-th structural and magnetic 
interfaces from the average height z n , respectively. 

To calculate the diffuse scattering from multilayers, in the spirit of Ref . Q we assume again 
for both E n (kj,/i) and E^(— k/, v) inside each layer the functional forms in the n-th layer, 
where z n < z < z n + Sz^ m \ n (x,y) for 5z( c>m y n (x, y) > 0, analytically continued to those in 
the (n + l)-th layer, where z n + Sz^ m ^ n (x,y) < z < z n for Sz( Ctm ) >n (x,y) < 0. Evaluating 
< kj, z/|T|kj, n > from Eq. (4.14) of paper 1^ without the first term on the right side and 
substituting Eqs. (|5.5|) and (|5.6|) give 

JV-l l,2(or cr,n) 3 

< k,, u\T\k h fi>=kij2 E E c; + \jj')F n (q^Uf)) , (5.7) 

n=0 jj' p=0 

where 





_ T (o) / 




v) T n+l,ji( k i,V), 




rp(0) f 
~ 1 n+l,j\ 


-k/, 






d(0) ( 


"k/ 


Z/ )^n+l,j'(k«, A 1 ), 


c 3 n+1 (jj") 


- r 


"k/ 


^i+ij'( k i,/i), 



(5.8) 



(Xn+l ~ Xn) E e *ja,n+l e j'a,n+l 



7;;: '(././') I V " ^ 

14 



x / / dxdye- iq ^ 1{jj ' )Sz ^ ix ' y) e- ici rP 



* / (2) _ (2) \ 

e ja,n+lv^o/3,n+l Xa/3,rJ e j'/3,n+l 

a/3 

x J J dxdye- iq ^ 1{jj ' )5zm -" ix ' y) e- i(l rp , (5.9) 

and 

C +1 0'iO = k fZjn+ i(j) - k iZtn+1 (f), q^Uf) = k fg>n+1 {j) - k r iZyn+1 (f), 

Substituting in Eq. (|4.1j) and carrying out the statistical average over the interfaces, we 
finally obtain 

/ rlrr \ diffuse £.4 N-l l,2(cr,7r) 3 

1 ] - E E E c; +1 uf)c;: +u (kk') 



dn 



167T 2 ^ 

n,n'=0 jj'kk' p,p'=0 



;^(c +i (jj / )))(^'*(<; +i (^ / ))>} ) 



(5.11) 



and 



where 



|Xn+i - Xn| 2 (E 



x (E e fc"^ / +l e fc , a.n'+l) Upc" 
a 

(Xaf3,n+1 ~ Xaj3,n) e j' P,n+Xj 



a/3 



e ka,n'+l{Xal3,n'+l Xaf3,ri ) e k' /3,n> +1 ) u mm 



a/3 



+ (Xn+1 - Xn)(E e ia,n+l e i'a,n+l 



^ (\~* I ( 2 )* ( 2 )* \ * \jjnn' , 

X VZ^ e ka,ri+l{Xa(3,n'+l ~ Xa/3,n' ) e k' f3,ri +1 J U cm + C.C. 

a/3 



= Ae 



X 



dXdYe-' 1 ^ 



■R ^^+ 1 (ii')<, 2 +1 *(^')^ i ',w(R) _ x ] . 



(5.12) 



(5.13) 



For circularly polarized incident x-rays, the scattering intensities without polarization anal- 
ysis for the outgoing beam can be also evaluated using Eqs. ()4.8|) and (|4.9j) . 
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We shall now illustrate numerical examples calculated for a multilayer with N interfaces. 
For the calculation of diffuse scattering intensities, we have assumed again the self-affine 
fractal interfaces for the height-height correlation functions.— For the cross-correlation func- 
tion between the ra-th and n'-th interfaces, the Schlomka et eZ.— expression was extended 
as 

l R l \ 2 V.n fW_\ 2h n'.n' !fn-*nH 



C HW (R) = ^±^(e + e v W j e HjT , (5.14) 

where C cc , C mm , and C cm are the cross- correlation functions for structural-structural, 
magnetic-magnetic, and structural-magnetic interfaces, respectively. Again, £u> represents 
the lateral correlation lengths, hu> the roughness exponents, and £± t ii> the vertical correlation 
lengths between different interfaces, respectively. 

Figure E] shows the results of calculations of x-ray resonant magnetic diffuse scattering 
intensities from a [Gd(51 A)/Fe(34 A)] 15 multilayer at the Gd L 2 -edge (7929 eV), which 
was used in Sec. IX of paper I.— Figure IHa)-(d) shows the dynamical calculations in the 
DWBA using Eq. ()5.11|) in Sec. V, whereas Fig. Hfe)-(h) show the kinematical ones 
using Eq. (j2.21|) in Sec. II. It can be clearly seen in the DWBA calculations that the 
anomalous scattering peaks^ in the rocking curves of (i+ — J_) intensities appear at the 
incident or exit angles corresponding to the positions of different order multilayer Bragg 
peaks, as shown in Fig. EJb)-(d). This has been observed experimentally by Nelson et al.fi 
but it has not been simulated theoretically because the kinematical calculation used by them 
cannot explain these Bragg-like peaks in the rocking curves even for the charge scattering 
intensities, as shown in Fig. 0fe)-(h). Nevertheless, kinematical calculations have been used 
widely because of their simplicity and good agreement of the overall features of the rocking 
curve with dynamical calculations. 

In order to compare our dynamical theory with experimental data, we used the same 
experimental data as that in Ref. Id, where the experimental data were fitted using the 
kinematical calculations. Figure El shows the measured sum (a) and difference (b) of opposite 
photon helicity rocking curve data (circles), as presented earlier in Fig. 4 of Ref. |9j, from 
a [Gd(53.2 A)/Fe(36.4 A)] i5 multilayer near the Gd L 3 -edge (7245 eV). The rocking curves 
were measured at the second (q z = 0.147 A -1 ) and the third (q z = 0.215 A -1 ) multilayer 
Bragg peaks. The lines represent the fits calculated in the DWBA using Eq. (I5.11|) . For the 
calculations, the charge and magnetic resonant scattering amplitudes near the Gd L 3 -edge 
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(7245 eV) were used as f c = 37.9 + 19. 8i and f m = —0.22 + 0.482, whose relationship to A 
and B defined in Eq. (3.3) of paper Ji was discussed in Sec. VIII of paper I.— Ferromagnetic 
layers were assumed to exist only near the Gd/Fe interfaces, and their layer thicknesses were 
7.8 A, as estimated in Ref. 0. From the best fit for both sum and difference intensities, 
we obtained the roughness amplitudes a c = 7.2 A and a m = 1.0 A, the lateral correlation 
lengths £ cc = 240 A and £ cm = 1000 A, the roughness exponents h cc = h cm = 0.3, and 
the vertical correlation lengths £_l )CC = 440 A and £_i_, C m = 670 A. When compared with 
the kinematical calculations presented as the solid lines in Fig. 4 of Ref. 0, the DWBA 
calculations in Fig. 5 show clearly that the anomalous scattering features indicated by the 
arrows in Fig. 5 can be explained well by the dynamical theory in the DWBA for both sum 
and difference intensities. 



VI. CONCLUSIONS 

We have shown in this paper and the preceding paper li that one can generalize the con- 
ventional theory of ordinary (Thomson) x-ray scattering from surfaces possessing self-affine 
structural roughness to the case of resonant magnetic x-ray scattering from surfaces or inter- 
faces of ferromagnetic materials possessing both structural and magnetic roughnesses. For 
this purpose, we have represented the deviations from a smooth magnetic interface in terms 
of "rough" magnetic interface, distinct from the structural interface (but possibly correlated 
strongly with it), with its own self-affine roughness parameters and parameters representing 
the correlation of the structural with the magnetic roughness height fluctuations. Compo- 
nents of the magnetization at the interface which are disordered on much shorter length 
scales are ignored in this treatment, as they will scatter at much lager q- values than those of 
interest here. The decrease of the in-plane averaged magnetization as a function of distance 
from the interface is taken into account by a form factor ip(q z ) which is the Fourier transform 
with respect to z of the derivative of graded average magnetization density, and a magnetic 
dead layer is taken into account by allowing for a possible separation A along the z-axis of 
the average structural and magnetic interfaces. In addition to magnetic roughness, magnetic 
domains can also give rise to offspecular scattering and their effect has also been included 
in the formalism. Formulae have been derived both in the Born approximation (BA) and 
the distorted-wave Born approximation (DWBA) for both single and multiple interfaces. 
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Numerical illustrations have been given for typical examples of each of these systems and 
compared with the experimental data from a Gd/Fe multilayer. 

We hope that the expressions given here and in the preceding paper will be useful 
in helping to analyze the rapidly increasing amount of magnetic x-ray scattering data cur- 
rently being accumulated from magnetic thin film and multilayer systems and in extracting 
meaningful parameters regarding to both the structural and magnetic roughness. This in- 
formation will help in the understanding of the magnetic and magnetotransport properties 
of these multilayered systems. The codes for the calculations in this paper and the preceding 
one are also available in C language by emailing to D.R.L. (drlee@aps.anl.gov). 
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APPENDIX A: EXPLICIT EXPRESSIONS FOR THE SCATTERING CROSS 
SECTIONS IN THE BORN APPROXIAMTION (BA) 

In order to obtain explicit expressions for the scattering intensities in Eqs. (|2.16|) or 
(I2.18J) . the polarization-dependent terms denoted by Eq. ()2.9|) should be calculated for a 
given scattering geometry— We here consider a common scattering geometry where k i; kf 
are both in the x — z plane (i.e., no out-of-plane scattering), as depicted in Fig. 1, and 

eM = &aiy) = y, 
e n (n) = xsin^j + zcos^, 

e n (u) = x(— sm Of) + zcosOf. (Al) 

/i, v represent the incident and final photon state, and 0{, Of are incident and scattered 
angles, respectively. 

Now we redefine G c and G m in Eq. ()2.9|) in terms of Gi^,3, which are more convenient 
for explicit calculations, by 

G c = (pi - P2 )G U G m = m m (BG 2 + tCG 3 ), (A2) 
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where 



Gi = (e*-e M ), 

G 2 = (e* v x e M ) • M, 

G 3 = (e:-M)(e M -M). 

Inserting Eq. ()A1|) in Eq. ()A3|) . we obtain 

= 1, G™ = G\ a = 0, GY = cos(9i + 0/), 



(A3) 



or 



0, = -M x cos % - M, sin 9 f} 



GY = M x cos 0i - M z sin U GY = M y sin(0; + 9 f ), 
Gl a = M 2 , GY = —M x M y sin 9 f + M y M z cos 9 f , 
GY = M x M y sin 9 { + M y M z cos 9 h 



GY = -M^ sm9i sin9 f + COS9.COS9 f + M Z M X sin^-fl/), 



(A4) 



where the first and second indices of the superscripts represent the polarizations of the inci- 
dent and final photon states, respectively. The offspecular scattering can be then expressed 
explicitly from Eq. ()2.18|) by 



da 

In 

da 
dtt 



2 

- 2 {\Pi ~ P2?S' CC + n 2 m \C\ 2 M*S' mm - 2( Pl - p 2 )n m Re[C*e-^ A }M 2 S' cm }, 



Q 



— J = ^n 2 n {\B\ 2 (M x cos9 f + M z sm9 f ) 2 + \C\ 2 M 2 (M x sm9 f - M z cos9 f f 



da s 

dn / 

da^ 

dn 



+ 2lm[BC*]{(M 2 x - M 2 )M y sm9 f cos9 f - M x M y M z cos(29 f )}}S' mm , 

2 

= ^n 2 m {\B\ 2 (M x cos 9i - M z 8^9^ + \C\ 2 M 2 (M X sin 9i + M z cos 9 { f 
+ 2lm[BC*]{(M 2 x - M 2 )M y sin 0< cos 9 t + M x M y M z cos(29i)}}S' mm , 
\ Pl - p 2 \ 2 cos 2 (9 t + 9 f )S' cc + n 2 m \\B\ 2 M 2 sm 2 (9, + 9 f ) 

+ |(7| 2 (M*sin0iSin0/ - M 2 cos 9 { cos 9 f - M z M x sm(9i - 9 f ) 
- 2\m[BC*}M y sm(9 l + 9 f ) 

x (Ml sin 6i sin 9 f - M 2 cos 9, cos 9 f - M Z M X sin(0 l - 9 f ) 
+ 2{ Pl - p 2 )n m [lm[B*e-^ A ]M y cos(9 i + 9 f ) sm(0, + f ) 
+ Re[C*e- iChA ] cos(9 t + Of) 

x (Ml sin 9i sin 9 f - M 2 cos 9 t cos 9 f - M Z M X sin(0; - 0/)) , v ,„ 



Si 



mm 



SI 



(A5) 
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12/2 2 \ * 16 

Replacing 5^- (i, j = c, m) by e" il3z ^~ z ^e~ ) anc l I| by -^r- 2 -, respectively, and setting 
9i = Of, the explicit expression of the specular reflectivity in Eqs. ()2.1fi|) and ()2.17|) can be 
immediately obtained. 

We also consider the case where the difference between the scattering intensities for right- 
(e + ) and left- (e_) circularly polarized incident x-rays is measured (mostly for ferromagnetic 
systems). From Eq. (jAl|) . the circular polarization vectors are given by 



e±0) = -^(^ ±«e 7r (/i)) = -^=(y ± i(xsin^ + zcos^)), 
e±(» = —f=(e a ± ie^u)) = -j=(y ±i(-±sin9 f +zcos0/)), 
and, inserting these in Eq. (jA3|) . we obtain 



(A6) 



Gt 



U 2 



Go 



GP 



Gt 



G— = i(l + cos(^ + 9f)), G+- = G- 1 + = l -{\- 008(9, + 9 f )), 



(G 2 
(G 2 
(G 



i(M x (cos9i + cos9 f ) + M z (sm9 f - sin^)) + M y sm(9i + #/)], 
i(M x (cos9i - cos^j) - M z (sin9 f + sin6> 4 )) - M y sin(6»; + 9 f )}, 



rr 



2 
1 

2 

-[i(M a ,Af y (sin0 i + sin0/) + M y M z (cos9i - cos9 f )) 



+ (-M 2 sin0j sin#/ + M 2 + M 2 cos^ cos% + M X M Z sin(^ - %))], 
1, 



(G 3 



[z(M E M 2/ (sin^ - sin0/) + M y M z (cos9i + cos0/)) 



+ (M 2 sin 0i sin % + M 2 - M 2 cos 0; cos 9 f - M X M Z sin(0j - 9 f ))] . 



(A7) 



Since the difference of the scattering intensities between positive and negative helicity of 
circularly incident polarization without polarization analysis for the outgoing beam can be 
evaluated as 



da" 



da" 

dn 



' da" 
dTL 



da 

In 



da" 

In 



da" 

dn 



(A8) 



inserting Eq. (A7) into Eqs. ()2.16|) and ()2.18|h we obtain 



' da\ (da 

dn),~ [dtt 



§{ (Pi - P2)n m [Re[ J B*e- i ^ A ](-l){M !r (cos^ + cos9 f cos(^ + 9 f )) 



+ M 2 (-sin 9i + sin 9 f cos(0i + 6 f ))\ 
+ Im[G*e- i ^ A ](-l){M :r M y (sin^ + sin^cos(^ + 9 f )) 



+ M y M z (cos9i — cos 9 f cos{9 i + 9 f ))} 



SL 
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+ n: 



^[BC^MxMl (cos Bi - sm9 f sin(^ + 9 f )) 
+ M Z M 2 (- sin 6i + cos 9 f sin(#, + 6 f )) 
- Ml sin 9i sin 8 f cos6 f + M z cos ^ sin % cos % 
+ Mr (cos 0< cos 2 0/ + sin 0/ sin(0; - 9 f )) 
+ .\ /;.U :: (cos Of sin(#j — 0/) — sin 0j sin 2 0/)}S'^ im >, 



(A9) 



and 



R. 



16ttV| 
? 2 4 



2(pi - p 2 )n m [Re[ J B*e- i ^ A ](-l)(M ;E cos 3 0, - M z sin 3 0,; 



+ 



+ n; 



Im[C'*e-^ A ](-l)(M E M 2/ sin0; cos 2 9 % + M y M z cos0; sin 2 0< 

Re^C*] [M X M 2 cos 0< cos(2^) + M 2 M 2 sin 9, cos(20<) - M 3 sin 2 0< cos 0* 



+ M 2 3 cos 2 0< sin 0; + M X M 2 Z cos 3 0< - M 2 M 2 sin 3 9 { 



(AlO) 



For a multilayer with N interfaces, we have assumed that resonant magnetic scattering 
amplitudes B n , C n of each layer have the same value B, C for all resonant layers, because 
these parameters depend primarily on each resonant atom itself. On the other hand, M n , 
n m ,n can have different directions and densities for each layer. We redefine again G Cjn and 
G m)n in Eq. (|2.20j) in terms of Gi,2,3 defined in Eq. flA3|) by 

G c ,n — (Pn+l ~ Pn)G\ = Ap n G\, 
G m ,n — iB [nm,n+lG2,n+l — ^771,11 ^?2,n] — C [n min+ iG 3tn+ \ — n mjn G 3jn ] 



BAG 2>n + iCAG 3 , n 



(All) 



and using Eq. (|A4|) the explit expressions for these terms in the cases of linearly polarized 
x-rays, i.e., aa, an, na, and irn scatterings, can be easily obtained. Since G±, AG 2j i, an d 
AG 3 ,i of all layers are real for linear polarizations as shown in Eq. (|A4|) . the offspecular 
scattering from a multilayer in Eq. (|2.21j) may be expressed in a more practical form as 

„2 N 



da 
dtt 



fj,— *v 



q 2 



^cc,nn 



\A Pn \*GiS' c 

+ {\B\*AG\ n + |Cf AG 2 n + 2lm[BC*]AG 2 , n AG 3 , n }s' mm , nn 
+ 2{{lm[A Pn B*]G 1 AG 2 , n - Re[Ap n C*]G 1 AG 3! n) cos(g 2 A n ) 

- (Re[A Pn B*]G 1 AG 2tn + Im[A / a n (7*]GiAG3,n) sin(g 2 A n )}^ 
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+ 2 {Re[Ap n Ap*,]G?cos(g z d n „,) + lm.[Ap n Ap* nl ]Glsm(q z d nnl )jS' CCtnn , 

n'>n 

+ {(\B\ 2 AG 2>n AG 2>n , + \C\ 2 AG 3in AG 3 , n> + lm[BC*]AG 2>n AG 3>n , 

- lm[CB*]AG 3in AG 2 ,n') cos(q z d nn >) 

+ (-Re[BC*]AG 2 , n AG 3 , nl + Re[C B*]AG 3jn AG 2jn ,) sm{q z d nn ,)}s' m 
+ {{(hn[Ap n B*]G x AG 2 , n , - Re[Ap n C , *]GiAG 3) „/) cos(g 2 (A„ + d nn ,)) 

- (Re[Ap n B*] G 1 AG 2 , n ' + MA^GxAGV) sin M A « + 



(A12) 

where {n <-> n'} refers to exchanging n and n' in the preceding term invloving with S' cm nn ,, 
d nn ' = z n — z' n , and S^ nn , = c, m) were defined in Eq. ()2.22|) . 

The difference of the scattering intensities from a multilayer between opposite helicities 
of circularly polarized x-rays can be also explicitly expressed from Eq. ()2.2H) and Eqs. 
(|S7 |) -lfgB |l as 



da 



da 



„2 N 



dn> I J .2 Yl 



dtt 



Q 2 



Re[BC*}AP{n,n-e % ,e f -Y )S' mn 

+ { {Re[Ap n B*}A?(n; 9 h 9 f ; T[) + lm[Ap n C*} A?(n; 9 U 9 f ; T' 2 )) cos(g 2 A„) 

+ (lm[Ap n B*]A?(n; 9 h 9 f ; T[) - Re[A Pn C*]A? (n; 9,, 9 f ; T' 2 )) sm(q z A n )}S' cr 



N r 

+ { (R4BC*] A[° } (n, n'; 6 U 9 f ; T ) + Re{CB*} A{? } (n, n'; 9 h 9 f ; T 5 )) cos(g,d nn 

n'>n 



+ (\B\ 2 AP (n, n'; i9 0,; T 3 ) + iCfA^ (n, n'; <? /; T 4 ) 
+ bn[SC*] Af°V, n'; i? 9 f ; T ) + lm[CB*]AP(n, n>; 9 h 9 f] T 5 )) sw{q z d nn ,)}S' mm>nn , 



+ 



U(Re[Ap n B*]AP(n'; 9 h 9 f ; r[) + Im[Ap n C*]AP (n'; 9 h 9 f ; T' 2 )) cos(q z (A' n + d nn >)) 
+ (lm[Ap„ J B*]AP(n / ;^^ /; r / 1 ) 

- Re[A Pn C*}A { r l \n'; 9 h 9 f ; T' 2 )) sm(q z (A' n + d nn ,))}S'„, 

+ [n n'}S' cm>n , n I , (A13) 



where 

Af, 0) (n, n'; 9 { , 9 f ; T p ) = n m , n+ in m . n ' +l T p (n + l,n'+ 1; 9 U 9 f ) + n m , n n m ^T p {n, n'; 9 i7 9 f ) 

r p (n+l,ri;6i,6 f ), 

A?{n-9 h 9 f -T' p ) = -n m ,n + iT' p {n + l-9 h 9 f )+n m , n T' p {n-9 h 9 f ), (A14) 
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and 



r o (n,n';0j,0/) = (M% cos6 f + M n z sin 9 f ) 

x ((Ml') 2 cos 0; cos 0/ + M£'M£' sm(6i - 6 f ) - (Ml') 2 sin 0; sin 0^ 
+ (M£ cos 0; - M 2 " sin 0;) (M y "') 2 

- M; sin(0i + 0/) (M X "'M;' sin 0/ - M;'M 2 "' cos 6 f ) , 

r' x (n; 0i,0 f ) = Ml cos 0; - M 2 n sin 0; + cos(0. 4 + 0/) (M£ cos 0/ + M 2 n sin 0/) , 
T' 2 (n; 9 h Of) = My sin 0, + M"M™ cos 0; + cos(0; + 0/) (M£M™ sin 0/ - M;M 2 n cos 0/) , 
r 3 (n,n';0<,0/) = M^ sin(0j + 0/)^ cos0/ + M 2 n sin0/) 

- M; sin((9i + 0/) (M;' cos 6 f + sin 6 f ) , 
r 4 (n, n'; 0,, 6 f ) = (M^) 2 (a#A/£ sin 0< + M;M 2 n cos 0;) 

- (M;) 2 (M^'My sin 0, + M;'M 2 n ' cos 0<) 

+ ((Af 2 n ') 2 cos 0^ cos 0/ + Ml' Ml' sin(0i - 0/) - (M™') 2 sin 0* sin 9 f ) 
x (m™M£ sin y - M^M'l cos 0/) 

- ({Ml) 2 cos 0j cos 0/ + Ml Ml sin(0i - 0/) - (Ml) 2 sin 0; sin 0/) 

x (Ml' Ml' sin 0/ - Ml' Ml' cos 0/) , 
r 5 {n,n';ei,e f ) = r {n',n;8i,9f). (A15) 

Again, replacing S' iinn , (i,j = c,m) by e ~* 9z(2i ' n ~^.™' ) e~ 5gz(CTl ' n+ ' 7 ^"' :) and % by 16 V° , re- 
spectively, and setting 0j = 0/, the explicit expression of the difference of the specular 
reflectivities with opposite helicities can be immediately obtained from Eq. (jAf 3(1 . 
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FIG. 1: Schematic of scattering geometry and sketch of the chemical (or structural) (z c (x,y)) and 
magnetic (z m (x, y)) interfaces, which can be separated from one another by an average amount A. 
Grazing angles of incidence (Si) and scattering (Of), the wave vectors kj and k/, and the photon 
polarization vectors of incidence (e M=(7j ^) and scattering (e u=rT;7T ) are illustrated. 
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FIG. 2: Schematic of rough structural and magnetic interfaces in a multilayer. 
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FIG. 3: Calculated x-ray resonant magnetic diffuse scattering intensities at the Gd L3-edge (7243 
eV) from Gd surfaces with different roughness parameters: (a) and (b) (a c , a m ) = (5 A, 3 A) and 
(£o Cm) = (800 A, 1500 A), (c) and (d) (a c , a m ) = (3 A, 5 A) with a 20-A-thick magnetically dead 
layer and (£ c , £ m ) = (1500 A, 800 A). Roughness exponent h = 0.8 was used for all structural 
and magnetic interfaces. In (c) and (d) the structural and magnetic interfaces separated by a 
magnetically dead layer were assumed to be completely correlated (circles), £j_ = oo. Solid (dashed) 
lines represent a — ► a (a — ► tt) scattering and open (filled) circles represent the positive (negative) 
values of the differences between I + and i_. Top panel: longitudinal diffuse scattering with an 
offset angle of 0.1°. Bottom panel: transverse diffuse scattering at q z = 0.2242 A -1 normalized 
to unity by the maximum diffuse scattering intensities to clarify the effect of lateral correlation 
lengths £ c and £ m . 
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FIG. 4: Calculated x-ray resonant magnetic diffuse scattering intensities from a [Gd(51 A)/Fe(34 
A)] 15 multilayer at the Gd L2-edge (7929 eV), which was used in Sec. IX of paper I»i The roughness 
amplitudes c C) Fe/Gd = 4.7 A, c Cj Gd/Fe = 3.6 A for charge interfaces and cr mFe / Gd = cr m G d / Fe = 4.2 
■^-j ^"m.Gd/Gd = 4.6 A for magnetic interfaces were used. For diffuse scattering, lateral correlation 
lengths (fee, £ cm , £ mm ) = (240 A, 1000 A, 1500 A), roughness exponents h cc = h cm = h mm = 
0.3, and vertical correlation lengths (£± lC e, £,±,cm, £±,mm) = (450 A, 670 A, 1400 A) were assumed. 
Top panel [(a)- (d)] shows the dynamical calculations in the DWBA using Eq. 1)5.11(1 in Sec. V, 
whereas bottom panel [(e)-(h)] shows the kinematical ones using Eq. (|2.21|) in Sec. II. (a) and (e) 
represent longitudinal diffuse scattering intensities with an offset angle of 0.1°, and (b)-(d) and (f)- 
(h) represent transverse diffuse scattering (rocking curve) intensities at the second- to fourth-order 
multilayer Bragg peaks, respectively. Thick (thin) solid lines represent a —* a (a —* ir) scatterings, 
and gray (black)-filled circles the positive (negative) values of the differences between I + and J_. 



The intensities of transverse scans are shown as a function of 
in order to illustrate the low q z scans better while q x = q z x 
radians. 
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FIG. 5: Measured sum [(J + + (a)] and difference [(/+ — I—), (b)] of opposite photon helicity 
rocking curve data (circles) at the second (q z = 0.147A -1 ) and the third (q z = 0.215A" 1 ) mul- 
tilayer Bragg peaks, which have been presented earlier in Fig. 4 of Ref. |9|- The lines represent 
the dynamical calculations in the DWBA using Eq. (|5.11|) and explain well the anomalous scat- 
tering features indicated by the arrows in both sum and difference intensities. The sample was a 
[Gd(53.2A)/Fe(36.4A)]i 5 multilayer, and the photon energy was tuned at 7245 eV (Gd L3-edge). 
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